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SUMMARY 


In  this  report,  all  the  data  pertinent  to  an  experimental  investigation  of 
the  effect  of  loading  system  rigidity  on  the  initial  buckling  load  of  un¬ 
reinforced  circular  cylindrical  shells  under  hydrostatic  pressure,  already 
reported  by  the  authors  (Reference  l)  in  condensed  form  are  presented. 
According  to  the  commonly  accepted  Tsien  criterion  (Reference  2),  the 
buckling  pressure  should  be  higher  in  a  rigid  system  than  in  a  soft  system. 
The  research  presented  in  this  report  denies  this  axiom  by  generating  strong 
evidence  that  the  stiffness  characteristic  of  the  loading  system  does  not 
alter  the  initial  buckling  stress,  which  fully  supports  the  results  obtained 
by  Kaplan  and  Fung  (Reference  3)  on  spherical  caps,  by  Krenzke  (Reference  4) 
on  plastic  spheres,  anu  by  Carlson,  Sendelbeck,  and  Hoff  (Reference  5)  on 
nickel-plated  spherical  shells. 

However,  in  the  case  discussed  in  this  report,  the  conclusion  is  based  upon 
a  statistical  interpretation  of  some  100  tests  on  very  nearly  identical 
vehicles ;  whereas  in  the  contemporary  researches,  a  more  limited  number  of 
vehicles  were  used.  Moreover,  it  concurs  with  the  similar  findings  by 
Horton,  Bailey,  Cox,  and  Smith  (Reference  6),  that  extensional  rigidity  of 
the  testing  machine  does  not  influence  the  initial  buckling  load  for  circular 
cylindrical  shells  loaded  in  axial  compression. 
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INTRODUCTION 


In  a  recent  paper  (Reference  I),  the  authors  reviewed  the  general  question 
of  the  influence  of  machine  stiffness  and  the  validity  of  the  Tsien  crite¬ 
rion  (Reference  2)  in  stiffness  computation.  The  present  report  presents 
all  the  test  data  for  one  series  of  experiments  discussed  in  that  presenta¬ 
tion.  The  question  herein  examined  is,  "Is  Tsien' s  criterion  pertinent  to 
the  instability  of  a  closed  thin -walled  vessel  under  external  pressure 
loading?"  Limited  experimental  data  by  Kaplan  and  Fung  (Reference  3)  and  by 
Krenzke  (Reference  4)  suggest  that  Tsien' s  criterion  is  not  applicable  in 
the  cases  of  shallow  spherical  caps  and  plastic  spheres .  The  present  study 
was  aimed  at  an  extensive  test  series  using  thin  cylindrical  shells  with 
end  caps  which  could  provide  the  basis  of  a  thorough  statistical  evaluation. 

The  results  obtained  were  in  full  agreement  with  the  conclusion  of  the 
researchers  referenced  above.  They  have  been  confirmed  further  by  subse¬ 
quent  researches  made  by  Carlson,  Sendelbeck  and  Hoff  on  spherical  shells 
Reference  (5). 
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EXPERIMENTAL  APPROACH 

Statistical  studies  of  the  behavior  of  structures  are  feasible  only  if  a 
large  number  of  nominally  identical  test  vehicles  are  available  at  reason¬ 
able  cost.  Thus,  the  vehicles  chosen  were  right -circular  cylinders  with 
end  caps.  Such  specimens  are  available  from  the  beverage  can  industry  and 
are,  as  previous  researches  have  demonstrated  (References  6  and  7),  remark¬ 
ably  consistent  in  character. 

Prior  experiments  had  illustrated  that  a  test  series  of  50  specimens  in 
each  of  two  environments  should  prove  to  be  adequate  for  statistical  dtudy; 
as  a  consequence,  this  number  of  vehicles  was  chosen.  The  choice  was 
further  reinforced  by  the  known  fact  that  cylindrical  shells  under  external 
pressure  loading  conditions  are  more  consistent  in  behavior  than  similar 
structures  under  axial  load.  This  is  clear  from  the  work  of  Strum  (Refer¬ 
ence  8)  and  from  the  extensive  study  of  Cleaver  (Reference  9)*  Thus,  the 
cylindrical  shell  appeared  to  be  a  potentially  better  choice  of  test  vehicle 
than  a  sphere,  since  it  is  well  recognized  that  minor  imperfections  in  such 
vehicles  have  extreme  significance  in  buckle  behavior.  A  great  quantity  of 
data  in  evidence  of  this  point  is  available  in  the  current  literature. 


DETAILS  OF  THE  TEST  SIECIMENS 


The  test  vehicles  were,  as  explained  in  the  preceding  section,  standard 
beverage  cans.  Such  cans  are  manufactured  in  a  fully  automated  process 
end  are  remarkably  consistent  with  regard  to  overall  geometry.  All 
specimens  were  made  from  the  same  batch  of  material.  Blanks  for  the  cylin¬ 
drical  bodies  were  punched  from  twice  cold-rolled  steel  sheet  and  then  were 
formed  into  a  cylindrical  shape  in  a  completely  automatic  process. 

The  overall  dimensions  were  as  follows: 

Diameter  -  2.63  inches 

Length  -  4.75  inches 

Thickness  -  0.0058  inch 

R/t  -  226 

In  the  interest  of  highest  consistency,  both  ends  were  capped  in  the  same 
automatic  machinery. 
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TEST  ENVIRONMENTS 


The  two  environments  chosen  were  obtained  by  change  of  fluid  within  the 
specimen.  In  one  case,  air  was  used;  in  the  other,  hydraulic  fluid.  In 
the  air  tests,  the  port  was  large;  in  the  hydraulic  fluid  tests,  small. 
Thus,  in  the  former  case,  the  volume  was  permitted  to  change  as  rapidly  as 
possible;  in  the  latter,  it  was  severely  restrained. 
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METHOD  OF  PRESSURE  APPLICATION 


For  ofce  soft  system,  the  method  of  application  of  load  was  as  delineated  in 
Figure  1.  It  is  clear  also  from  this  sketch  that  the  pressure  differential 
was  measured  by  means  of  a  standard  differential  manometer.  The  seal 
system  and  arrangement  of  the  specimen  in  the  test  apparatus  are  clearly 
apparent  from  Figures  2a  through  2  b.  Typically  specimens  collapsed  in  the 
same  manner  as  in  the  hard  system  (Figure  4). 

In  the  hard  system,  the  vessel  was  filled  with  hydraulic  fluid,  extreme 
care  being  exercised  to  ensure  that  there  was  no  air  entrainment .  The 
filled  vehicle  was  then  submerged  in  a  bath  of  hydraulic  fluid  to  mitigate 
the  effect  of  hydraulic  head.  This  system  was  then  connected  to  an  oil- 
mercury  differential  manometer.  The  full  details  of  the  pressure  loading 
system  are  clear  from  the  sketch  shown  in  Figure  3.  The  arrangement  of  the 
specimen  in  the  test  apparatus  is  portrayed  in  Figures  4a  through  4b.  A 
typical  buckled  specimen  is  depicted  in  Figure  4c. 


5 


■■ 

a.  Typical  Specimen  After  Collapse  From  External 
Pressure . 

Figure  2.  Continued. 
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Quick-Opening  Valve 


^TERMINATION  OF  INSTABILITY  PRESSURE 


The  determlne.tion  of  the  critic&l  pressure  was  readily  and  simply  made  by 
two  observers  who  continually  monitored  the  pressure  indicating  device. 
Collapse  in  either  system  was  readily  apparent  from  the  direction  change 
which  occurred  in  the  mercury  column.  However,  the  buckling  pressure  ob¬ 
tained  from  the  oil-mercury  differential  manometer  was  corrected  to  elimi¬ 
nate  the  effect  of  weight  of  oil  on  the  observed  pressure  head. 
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BEHAVIOR  UNDER  TESTING 


In  all  cases  the  failures  which  took  place  were  of  identical  character,  ini 
the  sense  that,  irrespective  of  the  stiffness  of  the  system,  the  shell 
distorted  frcm  a  circular  cross  section  to  a  square  cross  section.  However, 
in  the  soft  system  this  change  in  geometric  form  occurred  with  great  rapid¬ 
ity,  whereas  in  the  hard  system  it  was  very  slow.  The  final  buckle  pattern 
referred  to  above  is  apparent  from  Figure  2d  for  the  soft  system  and  from 
Figure  4c  for  the  hard  system. 
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DISCUSSION  CF  EXPERIMENTAL  DATA 

The  values  of  critical  pressure  for  the  soft  and  hard  systems  are  given  in 
Tables  £  and  II  in  the  appendix.  The  data  have  been  presented,  not  in  the 
order  of  tests,  but  according  to  magnitude  of  loads.  This  has  been  done  to 
simplify  the  analysis.  From  this  presentation,  it  is  simple  to  devise  cumu¬ 
lative  probability  distributions.  These  are  also  given  in  the  tables.  They 
have  been  graphically  portrayed  in  Figures  5  and  6.  It  is  Immediately  appar¬ 
ent  from  these  figures  that  both  distributions  are  essentially  normal.  As 
a  consequence,  the  data  are  analyzable  by  standard  statistical  procedures. 

The  mean  values  of  critical  pressure  are  derived  in  the  computations  given 
in  the  appendix.  It  is  shown  there  that  the  critical  levels  are  53.58  cm 
of  Hg  in  the  soft  system  and  53.73  cm  of  Hg  in  the  hard  system.  The  standard 
deviations  are  2.27  cm  of  Hg  and  2.11  cm  of  Hg,  respectively. 

From  the  close  identity  of  these  results,  it  might  be  stated  that  the  tests 
have  demonstrated  equality  of  performance  under  the  different  test  condi¬ 
tions.  However,  the  point  has  been  considered  in  depth  in  order  to  establish 
a  bond  of  confidence  within  the  confines  of  normal  statistical  testing. 

The  mean  buckling  pressures  have  been  compared  by  means  of  the  student's  "t" 
test.  This  comparison  is  made  in  the  appendix.  In  the  computation  given  in 
detail,  the  5$  level  of  significance  (  0!  a  0.05)  was  used,  and  the  standard 
derivations  were  assumed  to  be  equal.  The  calculations  show  that  for  a 
sample  size  of  50  and  a  5$  risk  of  accepting  a  false  hypothesis,  there  is  a 
95$  probability  of  detecting  a  difference  between  the  means  as  small  as 
1.6  cm,  or  about  3$  of  the  average  buckling  pressure.  The  t  test  statistic 
is  found  to  be  |  t|  =  0.349.  Thus,  since  0.349  <  t  Q  025.03  =  l-99»  we 
accept  the  hypothesis  that  there  are  no  differences  in 'the  mean  buckling 
pressures  obtained  from  the  variations  in  test  condition  at  the  5$  level  of 
significance. 
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conclusion 


The  results  presented  in  this  report  provide  substantial  experimental  evi¬ 
dence  that  the  stiffness  characteristic  of  the  loading  system  does  not 
affect  the  initial  buckling  stress  of  unreinforced  circular  cylindrical 
shells  tinder  hydrostatic  pressure.  This  conclusion  supports  the  experi - 
mental  findings  of  Kaplan  and  Fung;  Krenzke;  and  Carlson,  Sendelbeck,  and 
Hoff  on  the  buckling  of  sphere  caps  and  complete  spheres  by  external 
pressure  in  loading  systems  of  different  rigidities.  Thus,  the  commonly 
used  Tsien  energy  criterion  for  the  stability  studies  of  thin-walled  shells, 
which  predicts  a  higher  buckling  pressure  in  a  more  rigid  loading  system, 
must  be  an  improper  hypothesis  to  explain  the  observed  buckling  phenomenon 
of  shell  bodies. 
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APPENDIX 

STATISTICAL  ANALYSIS  OF  THE  EXPERIMENTAL  DATA 


In  this  appendix,  all  test  data  are  presented  and  analyzed.  The  individual 
critical  loads  for  the  tests  performed  in  the  soft  system  (air),  in  accor¬ 
dance  with  the  procedures  described  in  the  main  text,  are  listed  in  order 
of  magnitude  in  Table  I  and  graphically  portrayed  in  Figure  5 . 

The  buckling  loads  are  arranged  in  increasing  numerical  sequence.  In  order 
to  plot  all  data  values,  the  cumulative  probability  of  the  m  th  observation 
is  determined  from  m/(n  +  l) . 

The  data  pertaining  to  the  tests  made  in  the  hard  system  (oil)  are  likewise 
given  in  Table  n  .  This  information  is  graphically  portrayed  in  the  proba¬ 
bility  plot  of  Figure  6.  It  is  immediately  apparent  from  Figures  5  and  6 
that  the  distributions  are  essentially  normal;  therefore,  the  data  are 
analyzable  by  standard  statistical  methods. 

When  these  procedures  are  followed  with  the  data  relevant  to  the  soft  system, 

the  mean  buckling  load  ?  is  derived  as  follows: 

cr 


x 


2678.9 

50 


53, 578 -cm  Hg 


(1) 


2 

The  corresponding  sample  variance  (s  )  and  the  appropriate  standard  devia¬ 
tion  (s)  are  derived  from  the  equation 


2 

s 


(2) 


as  follows:  s  =  2875-751  -  2870.602 

and  s  =  2.268-cm  Jig 

while  the  coefficient  of  var  iation  v  , 
deviation  to  the  mean,  is  given  by 


=  5-149 

being  the 


ratio  of  the  standard 


v$ 


■4-  x  ioo$ 


2.268 

53-578 


x  100$  =  4.23$ 


(3) 
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TABLE  I.  CRITICAL  PRESSURE  FOR  50  NOMINALLY  IDENTICAL  CYLINDERS  TESTED 
UNDER  EXTERNAL  PRESSURE  -  AIR  PRESSURE  SYSTEM. 

nitude 


Load 

No. 

(n^) 

P 

cr 

Critical 
Pressure 
cm  Hg 

(*i) 

Cumulative 

Probability 

i 

m/(n  +  l) 

Load 

No. 

(m±) 

P 

cr 

Critical 
Pressure 
cm  Hg 

Cumulative 

Probability 

i 

m/(n  +  l) 

1 

48.6 

1.96 

26 

53.5 

51.0 

2 

49.1 

3.92 

27 

53.9 

52.9 

3 

49.4 

5.88 

28 

54.0 

54.9 

4 

50.3 

7.85 

29 

54.1 

56.9 

5 

50.6 

9.80 

30 

54.2 

58.9 

6 

50.7 

11.8 

31 

54.2 

60.8 

7 

50.9 

13.7 

32 

54.3 

61.6 

8 

50.9 

15-7 

33 

54.4 

64.6 

9 

51.1 

17.6 

34 

54.5 

66.6 

10 

51.5 

19.6 

35 

54.6 

68.5 

11 

51.5 

21.6 

36 

55-1 

70.5 

12 

51.7 

23-5 

37 

55.1 

72.5 

13 

52.1 

25-5 

38 

55.1 

74.5 

14 

52.2 

27.4 

39 

55.2 

76.5 

15 

52.5 

29.4 

40 

55-2 

78.5 

16 

52.6 

31.4 

4l 

55.5 

80.5 

17 

52.7 

33.3 

42 

55.8 

82.4 

18 

52.8 

35.1 

43 

55-8 

84.4 

19 

52.9 

37.2 

44 

56.2 

86.2 

20 

53.1 

39-2 

45 

56.6 

88.2 

21 

53.2 

4l.l 

46 

56.6 

90.1 

22 

53.3 

43.1 

47 

57.5 

92.1 

23 

53.3 

45.1 

48 

57.7 

94.0 

24 

53.4 

47.1 

49 

57-8 

96.0 

25 

53.5 

49.1 

50 

58.1 

98.0 

y p 

Lj  cr 

=  2,678.9 
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TABLE  II.  CRITICAL  PRESSURE  FOR  50  NOMINALLY  IDENTICAL  CYLINDERS 
TESTED  UNDER  EXTERNAL  PRESSURE  -  OIL  PRESSURE  SYSTEM. 
Values  arranged  in  order  of  magnitude 


Load 

No. 

(mi) 

Pcr 

Critical 
Pressure 
cm  Hg 

(xj 

Cumulative 

Probability 

* 

m/(n  +  l) 

Load 

No. 

(n^) 

P 

cr 

Critical 
Pressure 
cm  Hg 
(x1) 

Cumulative 

Probability 

% 

ra/(n  +  l) 

1 

49.1 

1.96 

26 

54.1 

51.0 

2 

49.3 

3.92 

27 

54.1 

52.9 

3 

49.3 

5.88 

28 

54.1 

54.9 

4 

51-.1 

7-85 

29 

54.1 

56.9 

5 

51-3 

9.80 

30 

54.4 

58.9 

6 

51.5 

11.80 

31 

54.6 

60.8 

7 

51.6 

13.7 

32 

54.6 

61.6 

8 

51.7 

15-7 

33 

54.7 

64.6 

9 

51.8 

17.6 

34 

54.7 

66.6 

10 

52.1 

19.6 

35 

55.0 

68.5 

11 

52.1 

21.6 

36 

55.1 

70.5 

12 

52.1 

23-5 

37 

55.2 

72.5 

13 

52.1 

25-5 

38 

55-4 

74.5 

14 

52.3 

27.4 

39 

55.4 

76.5 

15 

52.3 

29.4 

4o 

55-5 

78.5 

16 

52.3 

31.4 

4i 

55-6 

80.5 

17 

52.4 

33.3 

42 

55-9 

82.4 

18 

52.5 

35.1 

43 

55.9 

84.4 

19 

52.6 

37.2 

44 

55.9 

86.2 

20 

52.9 

39.2 

45 

56.2 

88.2 

21 

53.1 

4l.l 

46 

56.4 

90.1 

22 

53-1 

43.1 

47 

56.8 

92.1 

23 

53  A 

45.1 

48 

57.1 

94.0 

24 

53-8 

47.1 

49 

58.1 

96.O 

25 

53.9 

49.1 

50 

58.1 

98.0 

Y  p 

=  2,686.7 

L  cr 
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CUMULATIVE  PROBABILITY 


50 


51 


52 


53 


54 


55 


5* 


57 


BUCKLING  PRESSURE  t  Cm  Rg 


Figure  6.  Distribution  of  Buckling  Pressures  for  Thin  Circular 
Cylindrical  shells  When  Tested  in  a  Hard  System  -  Oil. 
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ThettKowtical.  curve  for  the  data  displayed  In  the  probability  plot  of 

X  =  P  +  (l/0f)y  (4) 

values  correspond  to  cumulative  probabilities  L~S  £  ES^'tST 

V^^^fto^*  X18vJichE?oSt,t°  6Stl“te  ^  Pam“tera 

P  *  x  -  53-573  cm.  Hg.  (5) 

'  °n  0.932  =  2*^  (6) 

mth  nS  50s  f'-TSS  ?‘“fard  deriation  ae  “  Of  sample  site, 

empirical  iL V  10,  Table  1.  2.  9,  p.  39).  Thua  the 

x  =  53*58  +  2.43y 
This  is  plotted  in  Figure  5. 

statistical8^ ormSa8 "  ^  eIr0ra  are  ^d'f?L  th/22 

a^xm^  — °  *  OOW  o(ym)  /n  (8) 

a  /n 

where  the  values  a(yj  /n  are  pure  numbers  independent  of  the  parameters 

(Reference  10  Tabled.  16,  p.  52).  ProbabiliiTvalues  ra^LT^  15 

r^°'5°sare^bta1^d  fr°”  syMKtrJr'  'rhe  control  curves  are  pitted  in’1 
Figure  5.  The  data  for  the  hard  system  (Oil)  is  treated  in  an  identical 
manner  and  the  various  parameters  of  importance  are  derived  as  follows: 


Buckling  -  p 
Load 


2686/ 


=  53*73^  cm.  Hg. 


Variance  =  s  =  2891.791  -  2887.343  =  4.448 
Standard  Deviation  =  s  =  2.110  cm.  Hg. 

Coefficient  of  variation  v$  =  100$  =  3.93$ 


The  coefficients  P  +  l/OC  are  given  by 

P  *  x  =  53.734  cm.  Hg. 

1/a  «  s/an  =  §^j§  =  2.264  cm.  Hg. 

The  line  defined  by  these  coefficients  is  plotted  in  Figure  6.  The  appropriate 
control  curves  are  derived  as  before.  The  necessary  values  of  o(x  )  being 
listed  below.  m 

a(xra)  =  ^  =  0.320  a(ym)/n  (ll) 

O/n 

These  control  curves  are  shown  in  Figure  6. 

It  is  readily  apparent  from  the  probability  plots  of  Figures  5  and  6  and 
their  appropriate  control  curves  that  both  sets  of  observations  can  be 
assumed  samples  from  normal  distributions.  The  standard  deviations  of  these 
distributions  are  extremely  close  and  therefore  the  means  ai-e  compared  using 
the  student  't'  test.  The  hypothesis  of  equality  is  examined. 

The  criterion  for  acceptance  of  this  premise  is 


(9) 

(10) 


I  t  I  s  V2.  n  +  n  -  2 
'  x  y 

where  a  is  the  level  of  significance 

and  n  +  n  are  the  sample  sizes, 
x  y 

Here  n  =  n  =50 
x  y 

n  =  n  -  2  =  98 
x  y 


(1£) 


(13) 


The  t  test  statistic  is  calculated  from  (Reference  11,  Table  7.2,  p.  171) . 


where  x  =  sample  data  from  soft  system 
y  =  sanple  data  from  hard  system 
x  =  53.578  y  =  53.734 

s2  =  5.1486  s2  =  4.448 
x  y 

nx  -  50i  ny  "  50 
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‘"(O'156’  dlsr  =°-3l« 

Choosing  a  5  percent  level  of  significance  a  table  of  percentage  points  of 
the  t  distribution  for  98  degrees  of  freedom  gives 

^0.025;  98  =  1,9S1* 


Since 


|  O.3I49  |  <  1.98^ 

the  hypothesis  of  equality  of  the  means  is  accepted  at  the  5  percent  level. 

An  indication  of  the  sensitivity  of  the  analysis  is  obtained  from  an  examina¬ 
tion  of  the  operating  characteristic  curve  at  this  level.  The  curve 
corresponding  to  a  sample  size  of  50  indicates  a  95  percent  probability  of 
detecting  a  difference  d  =  0.35*  (Reference  12,  Figure  6. 1C,  p.  129).  For 
this  test 

d  =  ^  1  (15) 

2a 

Using  the  average  sample  standard  deviation  from  the  two  tests  as  an  estimate 
for  a, 

|  p  -  p  |  =  0.70(2.189)  »  1.6  cm.  Hg.  (16) 

x  y 

Thus,  a  difference  between  the  means  as  small,  as  1.6  cm.  Hg.  or  approximately 
3  percent  of  the  average  buckling  load  could  be  detected  at  the  95  percent 
level. 
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